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ABSTRACT 


This report presents the computer program that has been 
used to analyze and design dual-throat electromagnetic dc con- 
duction pumps for the 5-kwe ZrH Reactor Thermoelectric System. 
In addition to a listing of the code and corresponding identification 
of symbols, a discussion of the bases for this analytical model is 
provided. 
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I. INTRODUCTION 


The 5-kwe Reactor Thermoelectric System (5-kwe System) is a compact, 
reflector-controlled, nuclear heat source, utilizing thermoelectric conversion 
to produce 5 kwe of electrical power continuously for 5 years. The prototype 
design of this SNAP system was completed by Atomics International, a division 
of Rockwell International. 

The 5-kwe System is illustrated in Figure 1; the principal components are a 
nuclear reactor, thermoelectric converter modules, dual-throat thermoelectric 
pump, heat rejection space radiator, volume accumulator units, interconnecting 
liquid-metal heat transfer primary and secondary coolant piping systems, and 
piping expansion joint units. 

Thermal energy produced in the nuclear reactor is transferred by liquid 
metal, NaK (binary eutectic 22 % sodium - 78% potassium alloy), circulated by 
the dual-throat thermoelectric pump through the reactor and thermoelectric 
modules of the primary loop, and through the thermoelectric modules and the 
radiator of the secondary loop. 

Previous SNAP system activities have included the design, fabrication, and 
testing of various types of dc electromagnetic Faraday- type pumps and pump 
systems. To facilitate the design of these pumps, the equations describing per- 
formance have been collected and set into computer codes. The ’'standard’ 1 
equations for predicting pump performance have been modified by experimental 
test results obtained from SNAP-type dc pumps. 

This report presents a computer model, established specifically for deter- 
mining minimum- weight dual- throat pump designs of the type shown in Figure 2 . 
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Figure 2. Prototype Thermoelectric Pump Assembly 







II. PUMP COMPUTER PROGRAM APPLICATION 


The computer program, printed in Section III, does not provide the optimum 
pump design, per se, but it furnishes the data required to make pump design 
point selections. DCPMP provides performance characteristics, and pump and 
pumping system weights, as a function of throat dimensions. Minimum pump 
weight and minimum pumping system weight do not occur at the same throat 
geometry, and neither geometry may provide the most desirable voltage, cur- 
rent, and efficiency characteristics. The pumping system weight includes the 
weight of the pump, the pump power supply (TE modules), associated support 
structure, and associated system components required to supply and reject heat. 
The weight of the thermoelectric power supply for pumps requiring from 1500 to 
1800 amp was set at 24 lb, and aweight assessment of 8 lb/kw of heat rejection 
was made for the associated heat supply and rejection equipment. The pump 
weight includes a penalty of 40% to provide the supporting structure for the pump 
and the TE modules. Selected pump performance parameters, such as voltage, 
current, power input, pump weight, and pumping system weight, are computed 
for many throat sizes. Plotting these parameters as a function of throat dimen- 
sions provides the throat and pump size at minimum pumping system weight, as 
well as selected key performance characteristics at this configuration. Selec- 
tion of a pump design point, using this technique, is discussed in Reference 1 
for the prototype pump. 

The computer program shown in this report makes a number of assumptions. 
These include: (1) identical current in each throat, (2) identical geometry in 

each throat, including transition pieces, (3) identical flow rate and pressure rise 
in each throat, and (4) identical flux density in each throat. In an actual pump 
design, the transitions may differ, the pressure - flow rate requirements may 
differ, and the flux density may differ. The requirements of each throat must be 
input separately, and the program run at the operating temperature of each 
throat, to determine which throat requires the greatest hydraulic pumping power, 
and therefore is dictating the pump size. In the prototype pump design of Refer- 
ence 1, the hydraulic pumping requirements were higher in the primary throat; 
thus, the primary throat determined the pump size for the 5-kwe Reactor TE 
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System. If the pumping requirements are significantly different for the two 
throats, it may be desirable to utilize different sized throats and magnet assem- 
blies; this would entail higher design and fabrication costs. 

The expressions used in this code to account for effects such as throat cur- 
rent leakage, interference effects between the current flow and the magnetic 
fields, and armature reaction effects are not universal expressions and are very 
dependent upon throat-magnet-bus geometry. These expressions were deter- 
mined from experimental data, and testing may be required to determine those 
parameters for pump designs that differ from that shown in Figure 2. 
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III. DCPMP AND PARAT COMPUTER CODES 


Following are the computer listings for the pump analytical model, DCPMP, 
and the temperature- sensitive parameters, PARAT. Table 1 identifies the sym- 
bols used in these programs. See Reference 6 for basic computer terminology. 


DCPMP 


100 
1 10 
1 20 
130 
1 40 
150 
1 60 
1 70 
180 

1 90 
200 
210 
220 
230 

2 40 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 


SCRATCH SPRITES 
SCRATCH s SEC ONL) : 

SCRATCH : TALLY : 

READ HO * B0,C8, C9 , D 

READ D 1 ,D2,D3,D4 

READ T,T1,T2,L3,Y 

DATA 645,1 1250, .02, .30*2.5 

DATA .265, .284, .321 , .290 

DATA .028, .022, .25,32, .025 

LET PI =1 .09 

INPUT: ONE :D5,U,R1 ,R2,R3,Z 
G0 T0 230 

INPUT*TW0:D5,U,R1 ,R2,R3,Z 
G0 T0 240 

F0R Q = 4 • 83 T 0 4.83 
F0R B= . 5 T 0 .7 STEP .1 
FOR A= 1 • 1 T0 1.5 STEP .2 
FOR X=2 T0 2.6 STEP .3 
IF A < 1 .5 THEN 300 
G0 T 0 320 
LET K =0 
G0 T 0 410 

IF A<= 1 .5 THEN 340 
G0 T 0 360 
LET K = 1 
G0 T 0 410 

IF A<=2 .25 THEN 380 

G0 T0 400 

LET K = 2 

G0 T 0 410 

LET K=3 

LET I1=R1*(A-K*T1 > / (B*X ) 

LET I2=R2*(A+2*T)/(2*T*X) 

LET I3=R2*(2*T+K*T1 )/(B*X> 

LET 14=2 • 6*R 1 /B 
LET 15=12*14/(12+14) 

LET Ml =5.714E“7*l5/( B* (15+11 )) 

LET M2 = 3 . 687E- 1 4/ ( D5*8 t2*(l5 + 1 1 )) 
LET M3 = . 96E- 1 4/ (D5*B*R 1 ) 

LET S=4*A*B/(3. 141 6*Dt2) 

LET N=8.64E4*Q/(U*(A+B) ) 
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510 IF N>2 • 3E3 THEN 540 
520 LET F 1 = 1 .1*1 6/N 
530 G0 T0 550 

540 LET FI =1 . 1*7 .9E-2*Nt-.25 

550 LET Cl =CC8+C9*C1 -S ) t2)*l .3E-3/(D5*(A-K*Tl ) t2*Bt2) 

560 LET C2=2*F1*<B+CA-K*T1 ))*(X+.4)*1 .3£-3/CDS*Bt3*(A-K*Tl )t3) 
570 LET M4=C 1 +C2 

580 LET I =2* ( (PI +M4*Q t2 ) * ( M2+M3 )*Q)t.5/Ml 
590 LET G=M1 *1 / C2* CM2+M3 ) *Q ) 

600 LET PI =M1 *G*I-(M2+M3)*GT2*Q-M4*Qt2 
610 LET El =6.452E-8*G*Q/(D5*B) 

620 LET VI =I5*(I*I 1 +E1 )/CI 1 +155 + 1*13 
630 IF Z=2 THEN 660 

640 PR I NTs PR I ME: A JBJX; I iGiQJPl JV1 iD5 

650 G0 T 0 710 

660 LET P2=P 1 

670 LET Q2=Q 

680 LET V2=V1 

690 LET D9=D5 

700 PRINT : SECOND :Q2 * P2 » V2* D9 

710 G0 T 0 720 

720 NEXT X 

730 NEXT A 

740 NEXT B 

750 NEXT Q 

760 IF Z=2 THEN 780 

770 G0 T0 220 

780 RESTORE: PRIME : 

790 REST0RE :SEC0ND: 

800 INPUT *PRIM£:A>B» X* I , G, Q, PI , V 1 , D5 

810 INPUT : SEC0ND : Q2* P2 * V2*D9 

820 LET V3=<V1 +V25/2 

830 LET V4=V 1 +V2+V3 

840 LET P7=I*V4 

850 LET P8=. 1 13*(P1 *Q/D5+P2*Q2/D9) 

860 LET 0=P8*1 00/P7 
870 LET LI =B+2*(T + Y 5 
880 LET F=5. 1 +L1+1 .23 
890 LET A1 =F*G*1 »05*A*X/B0 
900 LET L2=l .3*L1*G*1 .05/H0 
910 LET W1 =L2*A1*D1 
920 LET W2 =2 *T2*A*X*D2 
930 LET A3=R3*L3/CV3/I ) 

940 LET W3=L3*D3*A3 
950 LET D0=2*A*B/( A+B ) 

960 LET L=<D-DO 5/.21 

970 LET W4 = 2 *D4*T* (X*< A+B) +L* C A+B+l .571 *D) 5 +D4*K*T1*B*X 
980 LET C3= .03 l * ( A-K*T 1 )*B*(X+.4> 

990 LET C4=. 031*2. 0944+L* CD T2/4+0* (A+B 5 /6 .2832 +( A+B 5 »2/9 . 869) 
1000 LET W5=C3+C4 

1010 LET W6=D2*(L1 +L2+2*T2+5*A)*A1 *11 .25/19 
1020 LET W=1 .4*(W3+2*CW1 +W2+W4+W5+W6) ) 
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1030 LET W8=W+24 

1040 LET W9=<8*P7)/34 

1050 LET T8=W8+W9 

1060 PRINT A)B)X)I )G)W)T8)0J V4 

1070 PRINT : TALLY :A)BJXiWl;W2iW3iW4)W5lW6:W8JW9 
1080 G0 T0 800 
1090 END 


PARAT 


100 SCRATCH: 0NE 
1 10 SCRATCH : TW0 

120 PRINT “PRINT Y0UR PRIMARY AND SECONDARY THROAT TEMP IN DEG F M 
130 INPUT TM2 

140 LET R1 =(13.49+7. 463E-3*Tl+7. 205E-6*T1 t2)*l .OE-6 

150 LET R2=(27+.015*T1 >*1 .OE-6 

160 LET R3=(.58+.00148*(T1+T2)/2>*1 .OE-6 

170 LET Z=1 

180 LET D5=(54.27-.008349*T1 >/1728 
190 LET T = T1 +460 

200 LET U = 10t(. 6663 + 080. 26/T) - . 41 58*< CLG< T ) ) > 

210 PRINT : ONE: LNM( 10000) i D5 JU) R1 ) R2i R3 iZ) D 

230 LET Rl=(13. 49+7. 463E-3*T2+7.205E-6*T2t2)*l .OE-6 

240 LET R2=(27+.015*T2)*1 .OE-6 

250 LET D5=(54.27-.008349*T2)/1728 

260 LET T = T2 +460 

265 LET Z=2 

270 LET U=10t ( .6663+080. 26/T)-. 4158*CCLG(T) ) ) 

280 PRINT :TW0tLNM( 10000);D5|U)R1 )R2;R3iZ;D 
290 PRINT " " 

300 PRINT "DCPMP READY T0 RUN AT PRI AND SEC THROAT TEMP 0F" 

310 PRINT Tl>" & ", T2 

999 END 
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TABLE 1 


CODE SYMBOLS AND IDENTIFICATION 
(Sheet 1 of 3) 


Code 

Symbol 

Parameter 

Units 

Remarks* 

HO 

Maximum Field Intensity 

Oe 

Input 

BO 

Maximum Field Density 

G 

Input 

C8 

Throat Entrance Loss Coefficient 

- 

Input 

C9 

Throat Exit Loss Coefficient 

- 

Input 

D 

Piping Diameter 

in. 

Input 

D1 

Magnet Density 

lb/in. ^ 

Input 

D2 

Pole Piece Density 

lb/in. 3 

Input 

D3 

Bus Density 

lb/in. 3 

Input 

D4 

Throat and Splitter Density 

lb/in. ^ 

Input 

T 

Throat Wall Thickness 

in. 

Input 

T 1 

Splitter Thickness 

in. 

Input 

T2 

Pole Piece Length 

in. 

Input 

L3 

Bus Length 

• 

m. 

Input 

Y 

Magnet Air Gap 

in. 

Input 

D5 

NaK Density 

lb/in. ^ 

PARAT 

U 

NaK Viscosity 

lb/hr-ft 

PA RAT 

R1 

NaK Resistivity 

£i-in. 

PARAT 

R2 

Throat Wall Resistivity 

rZ-in. 

PARAT 

R3 

Bus Resistivity 

XZ-in. 

PARAT 

Z 

Throat Number (1 or Z) 

- 

PARAT 

PI 

Net Pump Pressure Head (Primary) 

! psi 

Input, PRIME 

P2 

Net Pump Pressure Head (Secondary) 

psi 

Input, PRIME 

Q 

NaK Flow Rate (Primary) 

lb/sec 

Input, PRIME 

Q1 

NaK Flow Rate (Secondary) 

lb/sec 

Input, PRIME 

B 

Throat Height 

in. 

Variable, PRIME 

A 

Throat Width 

in. 

Variable, PRIME 

X 

Throat Length 

in. 

Variable, PRIME 

K 

Number of Splitters 

- 

DCPMP 

11 

Resistance Through NaK 

a 

DCPMP 

12 

Resistance Through Side Walls 

\ a 

DCPMP 

13 

Resistance Through End Walls 
and Splitters 

a 

i 

DCPMP 
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TABLE 1 


CODE SYMBOLS AND IDENTIFICATION 
(Sheet 2 of 3) 


Code 

Symbol 

Parameter 

Units 

Remarks* 

14 

End Loss Leakage Resistance 
Through NaK 

a 

DCPMP 

15 

Total Bypass Leakage Resistance 

a 

DCPMP 

Ml 

See Section IV 

- 

DCPMP 

M2 

See Section IV 

- 

DCPMP 

M3 

See Section IV 

- 

DCPMP 

M4 

See Section IV 

- 

DCPMP 

S 

Area Ratio, Throat to Pipe 

- 

DCPMP 

N 

Reynolds Number 

- 

DCPMP 

FI 

Friction Factor for Throat 

- 

DCPMP 

Cl 

Entrance and Exit Hydraulic Losses 

psi 

DCPMP 

C2 

Throat Hydraulic Pressure Losses 

psi 

DCPMP 

I 

Total Current 

amps 

PRIME 

G 

Throat Flux Density 

G 

PRIME 

El 

Back emf 

V 

DCPMP 

VI 

Voltage Across Primary Throat 

V 

PRIME 

V2 

Voltage Across Secondary Throat 

V 

SECOND 

V3 

Voltage Drop Across Bus 

V 

Output 

V4 

Total Voltage Drop Across Pump 

V 

Output 

P7 

Power Input 

w 

Output 

P8 

Power Output (Hydraulic) 

w 

Output 


Efficiency 

% 

Output 

LI 

Distance Between Pole Faces 

in. 

DCPMP 

F 

Magnet Leakage Factor 

- 

DCPMP 

A1 

Magnet Cross-Sectional Area 

. 2 
in. 

DCPMP 

L2 

Magnet Length 

in. 

DCPMP 

W1 

Magnet Weight 

lb 

TALLY 

W2 

Pole Piece Weight 

lb 

TALLY 

A3 

Bus Cross-Sectional Area 

in. 

TALLY 

W 3 

Bus Weight 

lb 

TALLY 

DO 

Equivalent Throat Diameter 

in. 

DCPMP 
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TABLE 1 


CODE SYMBOLS AND IDENTIFICATION 
(Sheet 3 of 3) 


Code 

Symbol 

Parameter 

Units 

Remarks'" 

L 

Diffuser Length 

in. 

DCPMP 

W4 

Throat and Diffuser Weight 

lb 

TALLY 

C3 

NaK Weight in Throat at 70° F 

lb 

- 

C4 

NaK Weight in Diffusers at 70° F 

lb 

- 

W 5 

NaK Weight at 70° F 

lb 

TALLY 

W6 

Yoke W eight 

lb 

TALLY 

W 

Pump and Structure Weight 

lb 

Output 

W8 

Total Pump + TE Converter Weight 

lb 

TALLY 

W9 

Associated System Weight at 8 lb/kw 

lb 

TALLY 

T8 

Total Pump System Weight 

lb 

Output 


^Remarks 


Input 

PARAT 

DCPMP 

PRIME 

SECOND 

TALLY 

Output 


Read in via data statements 

Calculated by PARAT Program, and input upon request 

from Data Files ONE and TWO 

Calculated by DCPMP Program 

Primary Throat Characteristics Data File 

Secondary Throat Characteristics Data File 

Weights Data File 

Parameters selected for direct printout 
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IV. dc PUMP PERFORMANCE CHARACTERIZATION 


Performance characteristics of a dc electromagnetic-type pump can be 
determined from the equivalent electrical schematic and the standard pressure- 
flow equations of Reference 2. For the SNAP-type dc pump, this equivalent 
electrical schematic was modified as shown in Figure 3, by including the resist- 
ance of the splitter, the throat walls, and the bus that are electrically in series 
with the current flowing through the NaK in the throat. 



The symbols used in Figure 3 and in the following discussion are identical to 
those used in the computer program DCPMP (Section III), to facilitate cross 
checking between it and this discussion. All symbols are defined and listed in 
Table 1, in the order in which they appear in the program. Also, selected equa- 
tions are identified by line number (in parentheses), indicating their numerical 
location in the program. Therefore, from Figure 3, the voltage across the 
pump is: 


V1 = ITT 15 x ( IxI1 + E i) + i x 13 

where: 

TC _ 14 x 12 
15 = 14+12 


(620) 


(450) 
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From Reference 2, the back emf is: 


El = 6.452 x 10" 8 x G x Q/(D5 x B) 


(610) 


As a result of the pump program pump design, development, and testing, an 
additional term, to account for eddy current losses, was added to the pressure- 
flow relationship obtained from Reference 2. This term was derived in Refer- 
ence 3, and the value of the constant used with this expression was experimen- 
tally obtained. The complete pressure-flow relationship can now be expressed 
in the following form: 


PI = Ml x G x I - M2 x G 2 x Q - M3 x G 2 x Q - M4 x Q 2 


(600) 
(Equation 1) 


where: 


Ml x G x I = Gross pressure developed (psi) 

2 

M2 x G x Q = Pressure reduction due to back emf developed by fluid moving 
through a magnetic field (psi) 


M3 x G^ x Q 
M4 x 


= Pressure reduction due to eddy current (end) losses (psi) 

= Pressure reduction due to hydraulic losses in throat and in the 
entrance-exit transition sections (psi) 


Ml, M2, M3, and M4 are temperature-dependent constants, and are expressed 
as follows: 


Ml = (5.714 x 10- 7 /B) x I5 18 n 

(460) 

M2 = 3.687 x 10' 14 /[D5 x B 2 x (15 + 11)] 

(470) 

M3 = 0.96 x 1 0 — 1 4 /( D5 x B x Rl) 

(480) 

M4 = Cl + C2 

(570) 
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where: 


Cl = [C8 + C9 X (1 - S) 2 ] X 1.3 x 10 _3 /[D5 X (A - K x Tl) 2 X B 2 j (550) 

C2 = 2 x FI x [B + (A - K x T1)J x (X + 0.4) x 1.3 x 10" 3 / 

[D5 x B 3 x (A - K x Tl) 3 ] (560) 

The maximum magnetic flux density in the throat is obtained by differentiat- 
ing Equation 1 with respect to the flux, G. Thus 

G = Ml x I/[2 x (M2 + M3) x Q] (590) 

Substituting G back into Equation 1 and solving for I provides the current which 
will supply the required pressure and flow rate at this value of throat magnetic 
flux density: 


I = 2 x [(PI + M4 x O 2 ) x (M2 + M3) x Q] 1/2 /Ml (580) 

Higher flux densities (therefore longer magnets) will not provide greater pres- 
sure at this value of current and flow rate* 

The remaining equations in the program, for calculating power, efficiency, 
sizes, and weights, are of the standard textbook variety, and can be obtained 
from Reference 4. The voltage drop across the bus is fixed at one half of the 
total voltage across the two throats. Efficiency calculations include both throats 
and the total length of the current bus (^32 in. long). 

One of the key parameters in the analysis of the development pump (Refer- 
ence 1), in relating the experimental results with predicted results, is the cur- 
rent passing through the NaK beyond the ends of the active throat. This is ex- 
pressed as an end loss leakage resistance term in Line (440). 

14 = 2.6 x Rl/B 

The constant used in this expression was obtained from Reference 3, but the 
geometry of the experimental pump for which this constant was obtained was not 
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reported. The experimental data obtained from the development pump (Refer- 
ence 1) provided a much higher value (lower current loss) for this term, because 
the magnet length exceeded the bus length by an amount (2.75 to 2.3 in.) that pre- 
cluded much current from flowing outside of the magnetic field. Prediction of 
this term is one of the big uncertainties in the prediction of overall performance, 
and a careful examination of the magnet to bus relationship must be conducted 
by the pump designer. 

During the design effort on the NASA PLPA program (Reference 5), another 
term was added to the DCPMP program to account for armature reaction caused 
by distortion in the magnetic field intensity between poles due to electrical cur- 
rent flow. For uncompensated pumps, this effect can seriously impair pump 
performance at high currents (see Reference 2). This term was used in the com- 
puter program in the design of the prototype pump in Reference 1 until experi- 
mental data obtained from the development pump showed it was not required for 
a pump of this configuration and current level. It is presently not in DCPMP; but, 
for higher current pumps, may be considered. This term (K) was used in the 
following manner in the pressure-flow relationship: 


PI = Ml xGxIxK - (M2 + M3) x G 2 x Q - M4 x Q 2 


wh ere: 


K = 1 - ^ 277 ~^ 1 - Icoth B 1 


1 0 x G x B 


coth B1 - Tf 
Jt> 1 . 


B 1 = 


2 7T X Q X X 

10^ xAxBxRl x D5 
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V. PROGRAM USAGE AND OPERATION 


If this and all associated programs have been erased from the computer 
files, then the following sequence must be performed (if the programs are in the 
user's file, start with Step 3): 

1) Insert and save Programs DCPMP and PARAT, as listed in this report. 

2) Create dummy programs called ONE, TWO, PRIME, SECOND, and 
TALLY. DCPMP and PARAT will not run unless this is done. Just 
provide name and save each. 

3) Call PARAT and run. It will request primary and secondary throat 
temperatures. It computes NaK density and viscosity, and NaK, bus, 
and throat resistivities for each throat, storing the data in files ONE 
and TWO. The equations used in PARAT were supplied by NASA for 
the design effort conducted on the PLPA program (Reference 5). 

4) Call DCPMP and check all input values, Lines 140 to 270, and then 
run. The values stored with this program are for the prototype pump 
discussed in Reference 1 and shown in Figures 2 and 4. 

The code performs the following functions, in the order listed, for each 
input variable: 

1) Determines the number of splitters, Lines 280-400 

2) Computes resistances, Lines 410-450 

3) Computes temperature-dependent constants, Lines 460-570 

4) Computes current, Line 580 

5) Computes magnetic flux densities, Line 590 

6) Computes pressure (check of input pressure), Line 600 

7) Computes voltage across throat, Line 620. 

These equations are solved, and the results are stored in a data file, PRIME, 
for the primary throat, and then repeated for the secondary throat and stored in 
SECOND. 
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Additional equations are now solved in the following order for the primary- 
throat current and magnetic flux density: 

1) Bus voltage, Line 8Z0 

2) Total pump voltage, Line 830 

3) Input electrical power, Line 840 

4) Total hydraulic pumping power, Line 850 

5) Pump efficiency, Line 860 

6) Magnet dimensions, Lines 870-900 

7) Weights, Lines 910-1050. 

Selected parameters are printed out, Line 1060; for each input dimension and 
individual component weight data is stored in another data file, TALLY. 

After the throat dimensions are selected, these and the corresponding value 
of current and flux density are fixed in the program, Lines 580 and 590. The 
flow rate, Line 240, is varied through the range of interest to determine the 
head-flow relationship for the primary throat. The flux density is reduced ~20% 
to obtain the head flow relationship for the secondary throat. 
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